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Abstract 


This is the final report for a project funded by the Subsonic Assessment Program of the 
Atmospheric Effects of Aviation Project to implement an aerosol prognostic capability in the 
Global Modeling Initiative three-dimensional model. Though the original contract period was 
reduced from 3 years to 2 years and the funding was reduced by 50%, we have accomplished 
many of the original goals and are prepared to continue the work under other funding. The 
work involved collaboration between Atmospheric and Environmental Research, Inc. (AER) 
and the University of Michigan. The main accomplishment of the contract period was a thor- 
ough intercomparison of aerosol process algorithms between the AER microphysical code and 
the University of Michigan microphysical code. The AER code simulates nucleation, conden- 
sation, and coagulation of sulfate aerosols using a 40 bin size distribution. This approach is 
not practical for use in 3-D models because of its computational demands. The University of 
Michigan code calculates the mass loading and aerosol number concentration using a 2-mode 
distribution. A common set of algorithms which differ only in the method of resolving the 
aerosol size distribution has now been agreed to and is described in this report. Thus we are 
in a position to modify the 2-mode parameterizations to provide more realistic simulations. 
The AER two-dimensional aerosol model has been updated for consistency with the common 
algorithms and global calculations performed. 
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1 Project Overview 

The goal of this project was to develop an aerosol module to be used within the framework 
of the Global Modeling Initiative (GMI). The motivation came from the need to assess the 
atmospheric impact of sulfur compounds and particulate matter emitted by jet engines in the 
upper troposphere and lower stratosphere. The model development work was performed jointly 
by the University of Michigan and AER, using existing aerosol models at the two institutions 
as starting points. The GMI aerosol module was to be tested, evaluated against observations, 
and then applied to assessment of the effects of aircraft sulfur emissions as needed by the NASA 
Subsonic Assessment in 2001. The original scope of the work included the following tasks: 

1. Implementation of the sulfur cycle within GMI, including sources, sinks, and aqueous 
conversion of sulfur. Aerosol modules will be added as they are developed and the GMI 
schedule permits. 

2. Addition of aerosol types other than sulfate particles, including dust, soot, organic carbon, 
and black carbon. 

3. Development of new and more efficient parameterizations for treating sulfate aerosol nu- 
cleation, condensation, and coagulation among different particle sizes and types. 

4. Development of a 3-D database which can be used to evaluate results of the aerosol model. 

Note that the contract period was shortened from 3 to 2 years, with the total funding reduced 
to half of the proposed budget. Thus the performance objectives were modified over the original 
statement of work. Still, we were able to address most of the issues in the original proposal, 
though item (2) above was eliminated from the project scope. 

Sulfur source gas production and removal rates were delivered to GMI by Joyce Penner to 
begin the process of implementing the sulfur cycle in the 3-D model (Task 1). The GMI team, 
however, has been too busy to implement this yet. To accomplish Task 3, we performed inter- 
comparisons between the AER and University of Michigan aerosol aerosol models, comparing 
nucleation, condensation, and coagulation codes. This resulted in updates of both model codes, 
as common algorithms were agreed to. The final versions of the these algorithms are presented 
in this final report. We are now confident that the model differences are due only to the method 
of resolving the size distribution. Results of the model intercomparison exercises were presented 
at GMI meetings in January 2000 and June 2000. A 3-D database suitable for model evaluation 
(Task 4) exists at the University of Michigan and is available for model comparisons when GMI 
implements our aerosol modules and generates aerosol distributions. 


2 Scientific Results 

2.1 Approach 

The AER 2-D sulfate aerosol model is described in Weisenstein et al. [1997]. It uses 40 
size bins by volume doubling from 0.39 nm to 3.2 /rm. Each size bin is transported separately, 
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assuming a single constant radius. Condensation rates are calculated for each bin size, and 
coagulation deals with the interaction of each bin size with all others. This model has been 
used for the calculation of supersonic aircraft perturbations to stratospheric aerosol surface area 
density {IPCC [1999], Kawa et al. [1999]) and for the study of volcanic perturbations to the 
stratospheric aerosol layer {^Weisenstein et al. [1997]; Weisenstein and Yue [1999]). Nucleation 
previously followed the classical formulation of Zhao and Turco [1995] but has been replaced 
with the Kulmala et al. [1998] parameterization for this study. The diffusion coefficient in the 
condensation algorithm has been updated according to Jacobson [1999], and the Kelvin effect 
has been added. 

The University of Michigan 2-mode aerosol model follows Kveidenweis and Seinfeld [1988], 
Harrington and Kreidenweis [1998a], and Harrington and Kreidenweis [1998b]. The term 
“mode” refers to an aerosol distribution characterized by its mean radius, number concentra- 
tion, and a geometric standard deviation. The mass and number concentration of each mode 
are calculated prognostically in the model. The mean radius is derived from mass and number 
concentration; the standard deviation ag of the lognormal distribution is assumed. The size 
distributions of the two modes may overlap, though each mode is transported separately. Mode 
1 generally represents fine, submicron sized particles resulting from recent nucleation events. 
Mode 2 represents coarse particles which sediment rapidly. The prognostic equations for mass 
and number concentration include terms related to homogeneous nucleation of new sulfuric 
acid/water particles, condensation of H2SO4 on preexisting sulfate aerosol, and coagulation of 
sulfate aerosol particles. 

It was decided at the beginning of this project that the modal approach, currently used by 
the University of Michigan model, would be applied to the GMI model. The reason for this is the 
efficiency of the scheme, which makes it practical for three-dimensional modeling. Only aerosol 
mass and number in each mode need to be transported. The size-resolving approach used by 
AER requires that the number density in each of 40 bins be transported. Adding non-sulfate 
aerosol types would require an additional 40 bins for each type. The size-resolving approach does 
have the advantage of more accurately dealing with aerosol sizes and the interactions among 
particles. This would be particularly important in regions of nucleation and in regions such 
as aircraft wakes where aerosols are added to the atmosphere. To address this issue, we have 
performed model intercomparisons to analyze and improve the performance and accuracy of the 
modal scheme. We also have the option of increasing the number of modes for situations where 
more accuracy is required. 

Because the AER and University of Michigan models differ in dimensionality and transport 
parameterizations, a meaningful direct comparison was not possible. Therefore we chose to com- 
pare individual aerosol processes within a box model. Realistic initial conditions were chosen, 
some including an initial aerosol distribution and some initialized with only gas phase H2SO4. 
Analysis of these process model intercomparisons allowed us to determine differences between 
the AER and University of Michigan microphysical modeling approaches. Differences due to 
nucleation, condensation, and coagulation algorithms were resolved, so that the remaining dif- 
ferences are due only to the different methods of resolving the size distribution. The common 
algorithms agreed to are described below. 
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2.2 Microphysical Algorithms 

2.2.1 Nucleation 

The formulation of homogeneous nucleation follows the parameterization of Kulmala et al. 
[1998], with the nucleation rate J (particles cm“^ s~^) given by; 

J = exp(25.1289A,„,/-4890.8AT,„,^/T- 1743.3/T 

-2.2A796N,uifRH + 7643.4x„,/T - 1.97l2xai/RH) (1) 

with 


= ln(iV„/7V„,,) (2) 

, T - 273.15 
273.15 ■ 

Here is the ambient sulfuric acid vapor concentration, iVa ^ is the sulfuric acid vapor con- 
centration needed to produce a nucleation rate of 1 cm“^s~\ T is temperature (K), Xai is acid 
mole fraction of the embryo, and RH is fractional relative humidity. ,, is parameterized as: 



= exp(-14.5125 -I- 0.1335T - 10.546272// -|- 1958ARH/T). (4) 

The H2SO4 mole fraction in the critical nucleus is given by 
, 0.0154/2^1 

Xai = 1.2233 - — + 0.0102 In - 0.0415 In + 0.0016T (5) 

where RA denotes fractional relative acidity, and Nav and are sulfuric acid and water vapor 
concentration (cm“^), respectively. The embryo radius r* is obtained by solving 

. 2cru,„ 

0 ( 6 ) 

where 




I dp X a 
P dX a 




(7) 

( 8 ) 


after the composition (a:„/) is known. Here is the mass fraction of acid in the embryo, 
is Avogadro s number, ks is the Boltzmann constant, and p and o are the density and surface 
tension of the droplet. is the equilibrium concentration of water vapor over the surface 

of a solution with composition Xai- The parameterizations are valid over the temperature range 
from 233 K to 298 K and from 10 to 100% relative humidity. For use near the tropopause and 
in the lower stratosphere, we are awaiting an updated parameterization from Kulmala. 
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2.2.2 Condensation 


The condensation rate Gi onto Ni particles of radius n is given by 


( 9 ) 


with 


fi = 


1 + 


/ l.33j^ni + 0.71 4(1 -g) \ 

1 + Kui ^ 3tt / 


Kui = 


K 

Ti ' 


( 10 ) 

( 11 ) 


Ds is the diffusion coefficient for H 2 SOA molecules in air. AT, and iVf“‘ denotes the number 
density of H 2 SOi in the gas phase and the equilibrium H 2 SO 4 number density over the surface 
of sulfate aerosol particles, respectively. includes solution and curvature (Kelvin) effects 

[Jacobson, 1999]. /3 is the correction factor for non-continuum effects and imperfect surface 
accommodation. Two different values for the accommodation coefficient a are currently tested: 
g=l and a=0.04. in the equation for the Knudson number Kui denotes the mean free path 
for H 2 SO 4 molecules in air. 

In the two mode model, the particle volume mean radius, r„„/, times a correction factor is 
substituted for ri above 


Ti = ryoiexp{-ln‘^ag). 


( 12 ) 


This corrects for the effect of the size distribution on the condensation rate of the mode. However, 
the size dependence of the Knudson number Kn{ and of the supersaturation over a curved surface 
are ignored. Potential improvements to the two mode treatment of condensation could include 
a new parameterization of the correction factor or performing the condensation calculation for 
multiple bins within each mode. 


2.2.3 Coagulation 

We only consider Brownian coagulation. For the coagulation kernel Kij between a particles 
of radius ri and a particle of radius rj, we apply the interpolation formula of Fuchs [1964]: 


K, 




4 7r{Di + Dj) {ri + Tj) 

( 'Tx + rj 4 {Di -f D j) \ 

U, + r, + (p? + 5|)V2 (c2 cj)V2 {n + Tj)) 


(13) 


with Di = 


ksT 

6 m]ari 


0.87 


1 + ATna,i(1.249 + 0.42 exp(— - — )) 


= 2 


{Tj + Ix)^ - [r] + /?)^/^ 
3n h 


- 2 ri 


(14) 

(15) 

(16) 
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(17) 


li = 


4 A 

7T Ci 


T^rrii ) 


(18) 


A:b and T denote the Boltzmann constant and the temperature. A„ is the mean free path of air 
molecules, j]a the dynamic viscosity, and rrii the mass of a particle of radius Ti. 

In the 40 bin model, coagulation between two particles in the same bin results in one particle 
in the next larger bin. Coagulation between one particle and another particle from a smaller 
bin results in a fraction of the combined particle mass in the original larger bin and the next 
larger bin. In the two mode model, coagulation between two particles in the same mode results 
in a reduction in number density in that mode, but no mass change in that mode. Coagulation 
between particles in different modes results in mass transfer to the larger mode. A scheme to 
shift mass between modes as coagulation occurs has been implemented but needs further testing 
and refinement. 


2.3 Box Model Intercomparison 

The most important milestone in the first year of our project involved a model intercompar- 
ison of the two mode model and the 40 bin model. We first performed the comparisons using 
four sets of initial conditions, designated Box 1 through Box, 4. Relative humidities varied from 
0.66% to 95%, temperatures from 208 K to 263 K, and pressure levels were either 200 or 500 
mb. Each model used these initial conditions without transport and reported initial tendencies 
and time evolution of the aerosol distribution for a 24 hour period. Tables 1-4 show the initial 
conditions and initial rates for nucleation, condensation, and coagulation (Boxes 1 and 2 have 
no initial condensation or coagulation because there was no initial aerosol). This exercise re- 
sulted in many algorithm updates for both models as standard formulations were agreed to and 
differences resolved. We are now confident that the differences between model results derive 
solely from the different ways of treating the size distribution. 

The Figures 1-4 represent the time evolution of aerosol mass, number, nucleation rate, and 
condensation rate over 24 hours for each of the four initial conditions. Results from the two 
mode model (labeled UMICH) and the 40 bin model (labeled AER) are shown with values 
of O' in the condensation equation of 1.0 and 0.04. The value of a used in the condensation 
equation has generally been 1.0, which implies that all collisions between gas phase H2SO4 and 
aerosol particles result in particle growth. With a value of <y of 0.04, the condensation rate is 
dramatically reduced and the time required for gas-to-particle transfer is increased from 1 hour 
to 30 hours in Box 1 and from 1 hour to 12 hours in Box 4. Differences in the treatment of 
condensation (because of the different ways of representing the size distribution) are apparent 
and result in different number concentrations after integration. 

Figures 5-8 show the initial aerosol size distributions as a function of radius and those calcu- 
lated by each model during and at the end of the model integration for the case where oi—l.O. 
Box 1 results show that the two mode model accurately reproduces the size distribution in the 
situation where all aerosols are nucleated within the first minute and condensation ceases after 
2 minutes due to depletion of gas phase H2SO4. A single-mode log-normal distribution results 
after 1 hour due to coagulation, though the width of the distribution is somewhat narrower in 
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the 40 bin model than the two mode model. In Box 2, nucleation continues for 5 hours and 
the final aerosol mass is reached after 10 hours. In this case, the 40 bin model shows a skewed 
distribution after 1 and 5 hours due to continued particle generation in the smallest size bin, 
while the two mode model shows the same peak number density, but a symmetrical distribution. 

Box 3 was initialized with a three-mode size distribution, and subsequent changes in total 
aerosol mass were minor. Both models show no change in the size distribution over 24 hours, 
but the two mode model fits the initial distribution into a single mode, so the initial shape is 
not preserved. Box 4 had an initial aerosol distribution along with nucleation and condensation 
for the first hour, though the total aerosol mass increased by only 1% in the 24 hour period. 
The two mode model fit the initial aerosols into mode 2 and the newly nucleated aerosols into 
mode 1. The particle number densities in mode 1 after 1 and 5 hours are much less than in the 
40 bin model, and by 24 hours all mode 1 particles have been shifted into mode 2. 

We repeated the Box 4 calculations using a single lognormal distribution (rather than a 
superposition of three lognormals) as the initial condition. Figure 9 shows results of the Box 
4 calculation of aerosol mass, number concentration, nucleation rate, and condensation rate 
with the revised initialization. The AER results shown in red used the original initial aerosol 
distribution, which assumed 90% of the initial particles in a lognormal distribution with mode 
radius of 0.021 fim, 10% of the initial particles in a lognormal distribution with mode radius 
of 0.069 fim, and 0.04% of the initial particles in a lognormal distribution with mode radius 
of 0.351 fim. The green line represents a calculation performed with a single lognormal initial 
distribution with mode radius of 0.043 /xm. The UMICH model, since it fits the initial aerosol 
distribution into the larger of its two modes, can only represent the initial distribution as a 
single lognormal. Aerosol number concentration differs most between the models. AER model 
results using a single lognormal initial distribution (green lines) are considerably closer to the 
UMICH values than the standard AER results (shown in red). We find better agreement since 
the 2-mode model now accurately represents the initial distribution. However, compared to the 
AER model, the 2-mode model often overestimate condensation and underestimates coagulation 
rates. Some of this difference is likely due to the different timestepping schemes employed by 
the two models: an explicit scheme for AER and an implicit scheme for UMICH. 

An additional intercomparison exercise involved using a range of actual zonal mean and 
extreme conditions for January and performing 24 hour aerosol integrations. The calculations 
are for four altitudes (850, 500, 200, and 50 mb) and a range of latitudes. Figure 10 shows a 
scatter diagram of AER results vs. UMICH (labeled Bimodam) results as originally calculated 
in November 1999. Total particle number, nucleation over 24 hours, and condensation over 24 
hours are shown. The plots show that agreement between the models was reached only for a 
limited subset of conditions. Figure 11 shows current results after the models were updated to 
similar microphysics. Agreement is much improved and falls within reasonable limits for all but 
a few conditions. These results were presented at the International Conference on Nucleation 
and Atmospheric Aerosols by Michael Herzog. 

2.4 Enhancements to the AER 2-D Sulfate Aerosol Model 

Updates to the AER aerosol model have included implementation of new microphysical 
algorithms to match those agreed to for GMI and improvements to the 2-D transport and 
chemistry. The AER 2-D sulfate aerosol model and the AER 2-D chemistry transport model have 
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been combined into one model. Previously these were separate models, with the sulfate model 
using input OH and O3 from the chemistry model and the chemistry model using input surface 
area density from the sulfate model. With the models combined, perturbation calculations can 
be done with a single calculation rather than a series of two or more different calculations. For 
example, a Mt. Pinatubo calculation can calculate changes in sulfate surface area density which 
will impact ozone in the same calculation. While this currently has little eflfect on model results, 
it will allow for greater interactivity of sulfate and other trace gases in the future. For example, 
we could allow HNO3 to condense onto sulfate aerosols and generate PSCs of STS composition. 
This development work on the AER 2-D sulfate aerosol model was funded largely through our 
SAGE II contract (NAS 1-99096) but is of benefit to this contract also. 

Convection has been added to the AER 2-D model, using code and data files obtained from 
Victor Dvortsov \Dvortsov et al., 1998] and based on CCM3 convective fluxes. We found the 
convection to work well with the transport circulation of the GSFC 2-D model, as this transport 
is based on observational data. With the AER transport, which includes an enhanced 
term to partially account for convection, inclusion of the convective parameterization results in 
too much tracer transport between the surface and the upper troposphere. Figures 12 and 13 
show model-calculated aerosol extinction at 1.02 /xm and aerosol surface area density without 
(panel a) and with (panel b) convection. The calculated extinction with convection agrees well 
with SAGE observations for nonvolcanic periods. Aerosol surface area density is considerably 
enhanced by convection in the lower stratosphere, and may be somewhat too large compared to 
observations. 

The time stepping scheme in our 2-D aerosol model uses time splitting and performs nu- 
cleation before condensation. In cases with both large nucleation and condensation rates, the 
gas phase H2SO4 may be used up by nucleation and none left for condensation. Observational 
results indicate that under conditions with ample available surface area for condensation, the 
effects of nucleation are limited \Kulvnala et al., 2000]. We have explored this issue in the global 
2-D model by modifying the time stepping scheme. Figure 14a shows the annual average aerosol 
surface area density calculated by the model in its normal mode, where the nucleation step is 
performed before the condensation step. Figure 14b shows results with the condensation step 
performed first. Performing nucleation first results in 20-30% greater aerosol surface area density 
in the upper tropical troposphere. To help assess the accuracy of these schemes, we have subdi- 
vided the 1 hour time step into 40 substeps doing condensation and nucleation in sequence. The 
calculated surface area density with this scheme is shown in Figure 14c and falls between those 
shown in Figures 14a and 14b. This should represent a fairly accurate time-stepping scheme, 
though the computational time required is prohibitive. A new scheme has been devised which 
calculates the rate of change of gas phase H2SO4 due to nucleation and condensation (before 
modifying the aerosol number densities) and subdivides the time step if nucleation and/or con- 
densation will deplete all H2SO4 within a single time step. Results of this scheme are shown 
in Figure 14d. The computer time required for this scheme is quite modest because only grid 
points with large aerosol changes will perform a maximum of 50 substeps. 

The condensation and nucleation routines of the model have also been updated. Updates to 
the condensation rate include correction of the effective mean free path and diffusion coefficient 
according to Jacobson [1999], introduction of the Kelvin effect, and use of saturation vapor 
pressures from Ayers et al. [1980] and sulfuric acid activity coefficients from Taleb et al. [1996] 
as parameterized in Kulmala et al. [1998]. The impact of these changes on calculated annual 


8 



average surface area density for nonvolcanic conditions is shown in Figure 15, where 15a uses the 
old condensation routine and 15b the new routine. The most obvious impact is due to the Kelvin 
effect and results in a smaller surface area density at high latitudes and high altitudes. Use of 
updated vapor pressures in the condensation scheme has very little effect on calculated surface 
area density. It does have an impact, however, on nucleation rates. Figure 15c shows surface area 
density calculated with the new vapor pressures used in the Zhao and Turco [1995] nucleation 
scheme. Surface area density increases by 30-40% in the tropical upper troposphere with use of 
the updated vapor pressures. Figure 15d shows a calculation using the Kulmala et al. [1998] 
nucleation parameterization for grid points where the relative humidity is greater than 10% 
and the classical nucleation approach for all other cases. The Kulmala et al. parameterization 
includes the effects of hydration, which our classical approach does not. The nucleation rate 
in the troposphere is smaller with the Kulmala et al. parameterization, and thus surface area 
densities are smaller. 

2.5 Future Directions 

Implementation of a prognostic aerosol scheme for GMI should continue to be a high priority. 
The GMI model in the future will have the capability to adequately handle problems involv- 
ing both stratospheric and tropospheric chemistry and transport. Many problems involving 
atmospheric sulfate and aerosols (volcanic effects, aircraft impacts) will require both a realistic 
troposphere and a realistic stratosphere. If GMI funding continues, we are prepared to imple- 
ment the 2-mode aerosol module within GMI. The algorithms described here are the appropriate 
first step, but can be improved by increasing the number of modes (after model analysis to de- 
termine the optimum number), adding other aerosol types (dust, soot, black carbon, organic 
carbon), and developing pararaeterizations for the interaction of these aerosol types. 


3 Deliverables 

Source emission rates and chemical rate constants required to calculated gas phase sulfur 
chemistry were delivered to GMI by Joyce Penner. An aerosol database which will be used to 
compare three-dimensional calculated aerosol distributions with observations is available at the 
University of Michigan. Aerosol microphysical parameterizations are available as fortran code 
and can be delivered to the GMI team at Lawrence Livermore National Laboratory on short 
notice. 

The following presentations were related to this contract: 

1. given by Joyce Penner at the January 2000 GMI Meeting in Washington, DC [Penner et 
al., 2000]; 

2. given by Debra Weisenstein at the June 2000 GMI and AEAP Meeting, Snowmass, CO 
[Weisenstein et al., 2000]; 

3. given by Michael Herzog at the International Conference on Nucleation and Atmospheric 
Aerosols, Rolla, MO [Herzog et al., 2000a; 2000b]. 
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Table 1; Initial Conditions and Initial Rates for Box 1 


Quantity 

Units 

UMich Value 

AER Value 


Initial Conditions 


Pressure 

[mb] 

500 

500 

Temperature 

[K] 

262.96 

262.96 

Rel. humidity 

[%] 

86.57 

86.57 

Rel. acidity 

[%] 

13.10 

13.11 

H 2 O gas 

[molec/cm^] 

6.719x10^® 

6.723x10*° 

H 2 SO 4 gas 

[molec/cm^] 

7.097x10°® 

7.100x10°® 

Aerosol mass 

[molec/cm^] 

0.0 

0.0 

H 2 SO 4 mass fraction 

0.2078 

0.2078 

Aerosol density 

[g/cm^] 

1.1632 

1.1632 


Nucleation 


Nucleation rate 

[molec/cm^s] 

8.705x10°^ 

8.722 xl0°4 

Density of nucleus 

[g/cm^] 

1.562 

1.562 

Radius of nucleus 

[^m] 

5.006x10-°'* 

5.001x10-°'* 

Size of nucleus 

[molec/particle] 

3.199 

3.192 
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Table 2; Initial Conditions and Initial Rates for Box 2 


Quantity 

Units 

UMich Value 

AER Value 


Initial Conditions 


Pressure 

[mb] 

200 

200 

Temperature 

[K] 

233.93 

233.93 

Rel. humidity 

[%] 

0.66 

0.66 

Rel. acidity 

[%1 

27.89 

27.89 

H 2 O gas 

[molec/cm^] 

4.197x101^ 

4.199x10*° 

H 2 SO 4 gas 

[molec/cm^j 

1.407x10®^ 

1.408x10°^ 

Aerosol mass 

[molec/cm^j 

0.0 

0.0 

H 2 SO 4 mass fraction 


0.7486 

0.7486 

Aerosol density 

[g/cm^] 

1.6992 

1.727 


Nucleation 



Nucleation rate 

[molec/cm^s] 

0.590 

0.591 

Density of nucleus 

[g/cm^] 

1.805 

1.806 

Radius of nucleus 

[^m] 

4.272x10“°'* 

4.267x10"°'* 

Size of nucleus 

[molec/particle] 

3.199 

2.946 
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Table 3: Initial Conditions and Initial Rates for Box 3 


Quantity 

Units 

UMich Value 

AER Value 


Initial Conditions 


Pressure 

[mb] 

500 

500 

Temperature 

[K] 

208.12 

208.12 

Rel. humidity 

[%] 

6.366 

6.366 

Rel. acidity 

[%] 

100 

100 

H 2 O gas 

[molec/cm^l 

2.206 xlO^^ 

2.206x1048 

H 2 SO 4 gas 

[molec/cm^l 

2.603x10^5 

2.603 xl005 

Aerosol mass 

[molec/cm^j 

3.850xl0«7 

3.851 xl007 

Aerosol number 

[cm“^] 

8.503 

8.259 

Aerosol avg. radius 

[/im] 

5.670x10-02 

5.701x10-02 

H 2 SO 4 mass fraction 


0.6168 

0.6168 

Aerosol density 

[g/cm^] 

1.541 

1.587 


Nucleation 


Nucleation rate 

[molec/cm^/s] 

6.102x10-08 

6.102x10-08 

Density of nucleus 

[g/cm^j 

1.7905 

1.766 

Radius of nucleus 

[/im] 

3.064x10-04 

3.061x10-04 

Size of nucleus 

[molec/particle] 

1.032 

1.030 


Condensation 


H 2 SO 4 saturation 

[molec/cm^j 

3.056x10-08 

3.051x10-08 

Mean free path 

[cm] 

3.328x10-05 

3.327x10-00 

H 2 SO 4 gas diff. coeff. 

[cm^/ s] 

9.119x10-02 

9.124x10-02 

Condensation rate 

[molec/cm^/s] 

7.1213 

3.659 


Coagulation 


Particle velocity 

[cm/s] 

7.823 

7.707 

Coag. Kernel 

[cm^] 

1.137x10-0 

1.301x10-09 

Number Change 

[cm“^ s~^] 

-8.224x10-8 

-9.121x10-08 
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Table 4: Initial Conditions and Initial Rates for Box 4 


Quantity Units UMich Value AER Value 

Initial Conditions 


Pressure 

[mb] 

200 

200 

Temperature 

[K] 

218.05 

218.05 

Rel. humidity 

[%1 

95.43 

95.43 

Rel. acidity 

[%] 

100 

100 

H2O gas 

[molec / enU] 

1.110x10^^ 

1.110x10^'^ 

H2SO4 gas 

[molec/cm^j 

2.293x10'^^ 

2.293x10''^ 

Aerosol mass 

[molec/cnd] 

2.484x10*’^ 

2.485x10"^ 

Aerosol number 

[cm~^] 

411.08 

402.99 

Aerosol avg. radius 

[fim] 

5.700 xl0“^^‘^ 

5.701 xlO^*’-^ 

H2SO4 mass fraction 


0.1163 

0.1163 

Aerosol density 

[g/end] 

Nucleation 

1.091 

1.113 

Nucleation rate 

[molec/cnd/s] 

3.688 

3.689 

Density of nucleus 

[g/end] 

1.650 

1.636 

Radius of nucleus 

[fim] 

3.609x10“^’^ 

3.606x10“^ 

Size of nucleus 

[molec/particle] 

Condensation 

1.346 

1.343 

H2SO4 saturation 

[molec/cnd] 

1.458x10--^^ 

1.457x10"“ 

Mean free path 

[cm] 

8.718x10"'’^ 

8.715x10"”^ 

H2SO4 gas diff. coeff. 

[end/s] 

0.244 

0.245 

Condensation rate 

[molec/cnd/s] 

Coagulation 

4432 

2431 

Particle velocity 

[cm/s] 

9.517 

9.419 

Coag. Kernel 

[end] 

2.122x10“^ 

2.408 xl0"‘“ 

Number Change 

[cm”^ s”^] 

-3.587x10'^ 

-3.678x10"'“ 
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Figure 1: Modt'l calcTilations of aerosol evolution for Box 1 initial conditions from th(> UMICH 
(black lines) and tlu' AER (red lines) modc'ls. The change in ai'rosol mass, aerosol number 
conc(mtration. nucleation rate, and coiKhmsation rate are shown over 12 hours. Solid liiu's 
ix'juesi'iit valus of o in the condensation equation of 1.0, while dashed lines rei)r(‘.sent ft— 0.04. 
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Figuro 2: Model calculations of aerosol evolution for Box 2 initial conditions from the UMICH 
(black lines) and tin' AER (red lines) models. The change in aerosol mass, aerosol numbc'r 
conc('ntration. nuch'ation rate, and conden.sation rate an* shown over 24 hours. Solid lim's 


represc'ut valus of o in the' comh'usation eejuation of 1.0. whih' dashed lines ri'present n -().04. 
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Figure' 3: Model calculations of aerosol evolution for Box 3 initial conditions from the UMICH 
(black lines) and the AER (red lines) models. The change in aerosol mass, aerosol number 
conce'iitration. nucleation rate, and condensation rate are shown over 24 hours. Solid lines 
repn'sent valus of a in the condensation equation of 1.0, while dashed lines represent a— 0.04. 
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FiguiT 4: Modrl calcnlatioiis of aerosol evolution for Box 4 initial eomlitions from tli(> I'MICH 
(black liiK's) and tin' AER (red lines) inod<'ls. The change in m'rosol mass, aerosol mmiln'r 
conc('iitration. imcleation rate', and condensation rate are shown over 12 hours. Solid lim's 
rei)res('ut valus of o in the coiuh'nsation efiuation of 1.0, while dashed lines reim'sc'iit o— 0.04. 
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Radius (jam) 


Figure 5: Model caleulatioiis of tlie evolution of the aerosol size distribution for Box 1 initial 
conditions from th(> AER (top panel) and the UMICH (bottom panel) mod<ds. Results are 
shown after 1 hour. 5 hours, and 24 hours. 







Fif^un' 6: Model calnilations of th(> evolution of tlu' aerosol size distribution for Box 2 initial 
conditions from tli<' AER (toj) panel) aiul the FMICH (bottom ))anel) models. R('sults ate 
shown after 1 hour, 5 liours. and 24 hours. 
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Figure 7: Model ealeulations of the evolution of the aerosol size distribution for Box 3 initial 
conditions from the AER (top panel) and the UMICH (bottom jmnel) models. Results an* 
shown after 1 hour, 5 hours, and 24 hours, along with the initial ai^rosol siz(> distribution. 









Fi}i;ur(' 8: Modrl ralnilatioiis of the evolution of the aerosol size distribution for Box 4 initial 
conditions from the AER (to]) ])anel) and the UMICH (bottom panel) models. Results are 
shown aft(‘r 1 hour. 5 hours, and 24 hours, along with the initial a<'rosol siz(‘ distribution. 
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I'igun' 9: Mod('l calnilatioiis of aerosol evolution for Box 4 initial eonditions from tlu' UMICH 
(black lines) and the AER (red lines) models. The change in aerosol mass, aerosol numlu'r 
conc('ntration, nucleation rate, and condensation rate are shown over G or 12 hours. Solid lines 
repres(‘iit valnes of o in tlu' condensation equation of 1.0, while dashed lines rej>res('iit O—0.04. 
Ib'sults from tin' .A.ER model using oidy on(* lognormal mode for the initial condition are shown 
in green. 


24 




Box Model Comparison: 

first try 


Total Particle Number 

after 24h, January, extreme values 



Bimodam [#/cm'] 


Nucleation Condensation 



Figun' 10: Scattc'r diagiaiii of total particle number, nucleation rat(% and condc'iisation rat(' after 
a 24 hour calculation with the AER and UMICH models. Results rej)resent tin* first comparison 
of the models in Novc'inber 1999. 
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Box Model Comparison Today 

after corrections and improvements 


Total Particle Number 


after 24h, January, extreme values 



UMICH [#/cm'] 

Nucleation Condensation 



Figure 11: Scat tor diagram of total particle number, nucleation rate, and condensation rat(' after 
a 24 hour calculation with the AER and UMICH models. Results represent a comparison of tlu' 
current models after corrections and improvement to synchronize the microi)hysical treatment. 
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Figure 12: Calculated 1.02 fim aerosol extinction (10“^ km-1) for April under nonvolcanic 
conditions. The GSFC transport is used within the AER 2-D sulfate model, and panel (b) also 
includes the convective parameterization of Dvortsov et al. [1998]. 
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Figure 13: Calculated annual average aerosol surface area density cm— 3) under nonvol- 

canic conditions. The GSFC transport is used within the AER 2-D sulfate model, and panel 
(b) also includes the convective parameterization of Dvortsov et al. [1998]. 
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Figure 14: Calculated annual average aerosol surface area (/fm^/cm^) from the AER 2-D model 
using GSFC transport parameters and convection. Panel (a) is with nucleation performed before 
condensation, panel (b) is with condensation performed first, panel (c) is with condensation and 
then nucleation performed consecutively in a substep loop 40 times per hourly advective step, 
and panel (d) is with the new time stepping scheme. 
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Figure 15: Calculated annual average aerosol surface area (/im^/cm^) from the AER 2-D model 
using GSFC transport parameters and convection. Panel (a) is with our normal condensation 
scheme, panel (b) is with the updated condensation scheme, panel (c) is with the nucleation 
scheme updated with new saturation vapor pressures, and panel (d) is with the nucleation 
parameterization of Kulmala et al. [1998] used for relative humidities greater than 10%. 
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Project Overview 


A joint project between AER and the University of Michi- 
gan is underway to provide aerosol predictive capability in 
the Global Modeling Initiative (GMI) 3-D model. This 
project leverages off of the 3-D tropospheric two moment 
aerosol model at the University of Michigan and the 2-D 
stratospheric arosol model with 40 size bins at AER. The 
GMI aerosol module will initially operate with two aerosol 
modes for efficiency, but may later use 4 or more modes. 
Aerosol mass density and number density in each of the 
modes will be transported parameters, with condensation, 
nucleation, and coagulation modeled explicitly. Parameter- 
izations will be developed to account for the aerosol size dis- 
tribution in the condensation and coagulation calculations, 
using the 40-bin model as a guide. Eventually aerosol types 
other than sulfate could be added to the model, such as 
organic carbon, black carbon, dust, and sea salt. Modeling 
the interaction of aerosols with cirrus clouds and providing 
climate impact predictability are long-term goals. 
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Microphysical Model Descriptions 

The basic concept of the UMICH Aerosol Model follows 
Kreidenweis and Seinfeld [1988], Harrington and Krei- 
denweis [1998a], and Harrington and Kreidenweis [1998b]. 
The term “mode” refers to an aerosol subpopulation char- 
acterized by its mean size which is distinguished from a 
separately treated subpopulation of a different mean size; 
both are present in the same volume sample. Each mode is 
described by its mass and number concentration. For the 
size distribution within each mode, a lognormal distribu- 
tion with constant geometric standard deviation is assumed. 
The prognostic equations for mass and number concentra- 
tion include terms related to homogeneous nucleation of 
new sulfuric acid/water particles, condensation of H2SO4 
on preexisting sulfate aerosol, and coagulation of sulfate 
aerosol particles. 

The AER 2-D sulfate aerosol model is described in Weisen- 
stein et ai [1997]. It used 40 size bins by volume dou- 
bling from 0.39 nm to 3.2 /im. Each size bin is transported 
separately, assuming a single constant radius. Condensa- 
tion rates are calculated for each bin size, and coagulation 


deals with the interaction of each bin size with all others. 
This model has been used for the calculation of supersonic 
aircraft perturbations to stratospheric aerosol surface area 
density (/P (7(7 [1998], Kawa et al. [1999]). 

Nucleation 

The formulation of homogeneous nucleation follows the pa- 
rameterization of Kulmala et al. [1998]: 

J = exK25.1289A^,^//-4890.8Af,^///T- 1743.3/T 

-2.2m5NsuifRH + 7643.4xa//T - I.m2xai/RH) 

with 

Nsulf = lnNa/Na,c 

, , T- 273.15 

273.15 

Na is the ambient sulfuric acid vapor concentration, Na^c 
is the sulfuric acid vapor concentration needed to produce 
a nucleation rate of 1 cm“^s“^, T is temperature (K), Xai 
is acid mole fraction of the embryo, and RH is relative 
humidity. ^ ^nd Xai are also parameterized 
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Condensation 


The condensation rate G^ onto Ni particles of radius is 
given by 

Gi = 4 7T /? Ds {Ns - ""0 Ti Ni 


with 




Km = 


1 + 

A. 


/L33Xn, + 0.71 , 4(1 - a) 

^ 1 + Kfii 


-1 


+ 


3a 


Km 


Ti 


Dg is the diffusion coefficient for H 2 SO 4 , molecules in air. 
Ns and denotes the number of H 2 SO/^ molecules in 

the gas phase and the equilibrium H 2 SO 4 molecules over 
the surface of sulfate aerosol particles, respectively. N^^^ 
includes solution and curvature (Kelvin) effects [Jacobson^ 
1999]. P is the correction factor for non-continuum effects 
and imperfect surface accommodation. Two different values 
for the accommodation coefficient a are currently tested: 
a=l and a=0.04. Ag in the equation for the Knudson num- 
ber Km denotes the mean free path for H 2 SO 4 molecules 
in air. 

In the two mode model, the particle volume mean radius, 


r^, times a correction factor is substituted for above 

Ti — r^exp{—lTi?ag) 

This corrects for the effect of the size distribution on the 
condensation rate of the mode. However, the size depen- 
dence of the Knudson number Krii and of the supersatura- 
tion over a curved surface are ignored. 


Coagulation 


We only consider Brownian coagulation. For the coagu- 
lation kernel we apply the interpolation formula of Fuchs 
[1964]: 




47t(A- + A) A + ^ 7 ) 


Ti + r,- 


4 (A + Dj) \ 

' (c?-f rf)V2(r- + r^) 


-1 




with Di 




keT 

67r?7ari [ 

Aa 


1 + -fCna,j(1.249 + 0.42exp(- 


0.87 




)) 


r,' 




3 Ti li 


1 

H — 

TTCi 
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ks and T denote the Boltzmann constant and the temper- 
ature. Aa is the mean free path of air molecules and rja the 
dynamic viscosity. 

In the two mode model, coagulation between two par- 
ticles in the same mode results in a reduction in number 
density in that mode, but no mass change in that mode. 
Coagulation between particles in different modes results in 
mass transfer to the larger mode. In the 40 bin model, co- 
agulation between two particles in the same bin results in 
one particle in the next larger bin. Coagulation between 
one particle and another particle from a smaller bin results 
in a fraction of the combined particle mass in the original 
larger bin and the next larger bin. 

Box Model Intercomparison 

The first step in our project involved a model intercompar- 
ison of the two mode model and the 40 bin model. Four 
sets of initial conditions were chosen and designated Box 1 
through Box 4. Each model used these initial conditions 
without transport and reported initial tendencies and time 
evolution of the aerosol distribution for a 24 hour period. 


Tables 1-4 show the initial conditions and initial rates for 
nucleation, condensation, and coagulation. This exercise 
resulted in many algorith updates for both models as stan- 
dard formulations were agreed to. The differences between 
model results now derive solely from the different ways of 
treating the size distribution. 

The figures represent the time evolution of aerosol mass, 
number, nucleation rate, and condensation rate over 24 
hours for each of the four initial conditions. Results from 
the two mode model (labeled UMICH) and the 40 bin model 
(labeled AER) are shown with values of a in the condensa- 
tion equation of 1.0 and 0.04. Differences in the treatment 
of condensation (because of the different ways of represent- 
ing the size distribution) are apparent. Also shown are the 
initial aerosol size distributions as a function of radius and 
those calculated by each model during and at the end of the 
model integration. 
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Table 1 ; Initial Conditions and Results for Box 1 


Quantity 

Units 

UMich Value 

AER Value 


Initial Conditions 


Pressure 

[mb] 

500 

500 

Temperature 

[K] 

262.96 

262.96 

Rel. humidity 

[%] 

86.57 

86.57 

Rel. acidity 

[%] 

13.10 

13.11 

H 2 O molecules 

[molec/cm^l 

6.719E+16 

6.723E+16 

H 2 SO 4 gas 

[molec/cm^l 

7.097E+08 

7.100E+08 

Aerosol mass 

[molec/cm^] 

0.0 

0.0 

H 2 SO 4 mass fraction 


0.2078 

0.2078 

Aerosol density 

[g/cm^l 

Nucleation 

1.1632 

1.1632 

Nucleation rate 

[molec/cm^/s] 

8.705E+04 

8.722E+04 

Density of nucleus 

[g/cm^l 

1.562 

1.562 

Radius of nucleus 

[/zm] 

5.006E-04 

5.001E-04 

Size of nucleus 

[molec/part] 

3.199 

3.192 
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Table 2: Initial Conditions and Results for Box 2 


Quantity 

Units 

UMich Value 

AER Value 


Initial Conditions 


Pressure 

[mb] 

200 

200 

Temperature 

[K] 

233.93 

233.93 

Rel. humidity 

[%] 

0.66 

0.66 

Rel. acidity 

[%] 

27.89 

27.89 

H 2 O molecules 

[molec/cm^j 

4.197E+13 

4.199E+13 

H 2 SO 4 gas 

[molec/cm^j 

1.407E+07 

1.408E+07 

Aerosol mass 

[molec/cm^j 

0.0 

0.0 

H 2 SO 4 mass fraction 


0.7486 

0.7486 

Aerosol density 

[g/cm^j 

1.6992 

1.727 


Nucleation 



Nucleation rate 

[molec/cm^s] 

0.590 

0.591 

Density of nucleus 

[g/cm^j 

1.805 

1.806 

Radius of nucleus 

[/im] 

4.272E-04 

4.267E-04 

Size of nucleus 

[molec/part] 

3.199 

2.946 



Table 3: Initial Conditions and Results for Box 3 


Quantity 

Units 

UMich Value 

AER Value 


Initial Conditions 


Pressure 

[mb] 

500 

500 

Temperature 

[K] 

208.12 

208.12 

Rel. humidity 

[%] 

6.366 

6.366 

Rel. acidity 

[%1 

100 

100 

H 2 O molecules 

[molec/cm^j 

2.206E+13 

2.206E+13 

H 2 SO 4 gas 

[molec/cm^] 

2.603E+05 

2.603E+05 

Aerosol mass 

[molec/cm^l 

3.850E+07 

3.851E+07 

Aerosol number 

[cm“^] 

8.503 

8.259 

Aerosol avg. radius 

[/im] 

5.670E-02 

5.701E-02 

H 2 SO 4 mass fraction 


0.6168 

0.6168 

Aerosol density 

[g/cm^j 

1.541 

1.587 


Nucleation 



Nucleation rate 

[molec/cm^/s] 

6.102E-08 

6.102E-08 

Density of nucleus 

[g/cm^j 

1.7905 

1.766 

Radius of nucleus 

[/im] 

3.064E-04 

3.061E-04 

Size of nucleus 

[molec/part] 

1.032 

1.030 


Condensation 


H 2 SO 4 saturation 

[molec/cm^] 

3.056E-03 

3.051E-03 

Mean free path 

[cm] 

3.328E-05 

3.327E-06 

H 2 SO 4 gas diff. coeff. 

[cm^/s] 

9.119E-02 

9.124E-02 

Condensation rate 

[molec/cm^s] 

7.1213 

3.659 


Coagulation 



Particle velocity 

[cm/s] 

7.823 

7.707 

Coag. Kernel 

[cm^] 

1.137E-9 

1.301E-09 

Number Change 

[cm~^ s“^] 

-8.224E-8 

-9.121E-08 
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Table 4; Initial Conditions and Results for Box 4 


Quantity 

Units 

UMich Value 

AER Value 


Initial Conditions 


Pressure 

[mb] 

200 

200 

Temperature 

[K] 

218.05 

218.05 

Rel. humidity 

[%] 

95.43 

95.43 

Rel. acidity 

[%1 

100 

100 

H 2 O molecules 

[molec / cm^j 

l.llOE+15 

l.llOE+15 

HaSOj gas 

[molec/cm^l 

2.293E+06 

2.293E+06 

Aerosol mass 

[molec/cnd] 

2.484E+08 

2.485E+08 

Aerosol number 

[cm“'^] 

411.08 

402.99 

Aerosol avg. radius 

[^m] 

5.700E-02 

5.701E-02 

H 2 SO 4 mass fraction 


0.1163 

0.1163 

Aerosol density 

[g/cnd] 

1.091 

1.113 


Nucleation 



Nucleation rate 

[molec / cnd/s] 

3.688 

3.689 

Density of nucleus 

[g/cm^] 

1.650 

1.636 

Radius of nucleus 

[/mi] 

3.609E-04 

3.606E-4 

Size of nucleus 

[molec/part] 

1.346 

1.343 


Condensation 


H 2 SO 4 saturation 

[molec/cnd] 

1.458E-11 

1.457E-11 

Mean free path 

[cm] 

8.718E-06 

8.715E-06 

H 2 SO 4 gas diff. coeff. 

[cnT/s] 

0.244 

0.245 

Condensation rate 

[molec/cm^/s] 

4432 

2431 


Coagulation 


Particle velocity 

[cm/s] 

9.517 

9.419 

Coag. Kernel 

[cm^] 

2.122E-9 

2.408E-09 

Number Change 

[cm"^ s"^] 

-3.587E-4 

-3.678E-04 
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[10*^ molec/cmVs] molec/cm ] 


Box 1 no gas source 

Aerosol Mass 



time [min] 


Nucleation Rate (10s averages) 



time [min] 






[molec/cmVs] [lO" molec/cm'] 










[10^ molec/cm '] 





















Concentration (#/cm' ) Concentration (#/cm' ) 


Box 2 Size Distribution Evolution 



Radius (urn) 



Radius (}im) 
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I I I I ml i I i 1 1 ml i III mil I « « « ■ ml 








Concentration (#/cm ) Concentration (#/cm^) 


Box 3 Size Distribution Evolution 
















